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We investigate spin relaxation in a silicon double quantum dot via leakage current through Pauli 
blockade as a function of interdot detuning and magnetic field. A dip in leakage current as a function 
of magnetic field on a ~ 40 mT field scale is attributed to spin-orbit mediated spin relaxation. On 
a larger (~ 400 mT) field scale, a peak in leakage current is seen in some, but not all, Pauli-blocked 
transitions, and is attributed to spin-flip cotunneling. Both dip and peak structure show good 
agreement between theory and experiment. 



Electron spins confined in semiconductor quantum 
dots (QDs) are attractive candidates for quantum infor- 
mation processing [1 . Coherent manipulation of indi- 
vidual and coupled electron spin states has been mainly 
investigated in GaAs-based double QD (DQD) devices 
[2J-0] . However, nuclear spins of the host material cause 
decoherence of the electron spin via strong hyperfine cou- 
pling [5 . To reduce this effect, group IV materials, such 
as carbon, silicon (Si), and silicon-germanium (SiGe), 
have been investigated [6HTD] because their most abun- 
dant isotopes have zero nuclear spin. Silicon systems, in 
particular, have an advantage for future integration be- 
cause of their compatibility with conventional Si metal- 
oxide-semiconductor devices. 

Toward spin qubits in Si systems, it is necessary to 
understand the spin relaxation mechanism. Pauli spin 
blockade (PSB) pTJ[12] is a valuable tool for investigating 
spin relaxation in confined systems. In DQDs of several 
materials, the spin relaxation mechanism has been char- 
acterized by analyzing the leakage current in the PSB 
regime [13-16 , where hyperfine interaction and/or spin- 
orbit interaction dominate the spin relaxation. For Si 
systems, a PSB has been reported for a DQD in metal- 
semiconductor-oxide structures and an electrostatically 
formed DQD in Si/SiGe heterostructures [17j[T8]. How- 
ever, the relaxation mechanism in Si DQDs has not yet 
been experimentally clarified. More recently, magnetic 
field dependences of the leakage current in a PSB regime 
have been demonstrated in a pure Si DQD [19 , where a 
current peak was explained by field-dependent cotunnel- 
ing. 

In this Letter, we investigate leakage current in a 
PSB regime using a lithographically defined Si DQD. By 
changing magnetic field, we observed a dip of the leak- 
age current at zero magnetic field, presumably the re- 
sult of spin-orbit-mediated spin relaxation. In addition, 



magnetic field dependences at a different charge triple 
point exhibit a leakage current peak at zero magnetic 
field. This peak can be understood as a signature of 
spin-flip cotunneling processes. 

Figure [TJa) shows a schematic of a Si DQD. Three 
constrictions between the source (S) and drain (D), and 
five side gates were patterned by electron beam lithogra- 
phy on a 60-nm-thick (100) Si-on-insulator (SOI) layer, 
where the thickness of the buried oxide was 400 nm. Re- 
active ion etching was used to transfer the resist pattern 
onto the SOI, followed by formation of the gate oxide via 
thermal oxidation for 30 min at 1000 °C and low-pressure 
chemical vapor deposition (LPCVD). Then, a wide poly- 
Si top gate (TG) formed by LPCVD was used as an ion 
implantation mask for the formation of the n-type S and 
D regions. Finally, 300-nm-thick aluminum contact pads 
were formed by electron beam evaporation. Figure [TJb) 
shows a scanning electron microscope image of the de- 
vice, where the DQD is defined by tunnel barriers at the 
three constricted regions [20] . 

Electrons were attracted to the Si (100) surface by ap- 
plying a positive TG voltage, Vtg- Electrochemical po- 
tentials of the left and right QDs were modulated by 
applying voltages Vl and Vr to side gates L and R. The 
tunnel coupling between the two QDs was controlled by 
voltage Vc applied to side gate C. All measurements were 
carried out in a 3 He refrigerator with a base temperature 
of 250 mK. 

The honeycomb charge stability [Fig. [TJc)] reflects the 
formation of a DQD [21 j. Charging energies of the left 
and right QDs were estimated to be 10.7 and 11.0 meV, 
respectively, from the spacings of the Coulomb peaks, im- 
plying that the QDs have almost the same size. In addi- 
tion, from the distribution of the current peaks due to res- 
onant tunneling at triple point A in Fig. [TJc), the quan- 
tum level spacing, AE, of the left and right QD was esti- 
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FIG. 1: (color online) (a) Schematic of the silicon double 
quantum dot (Si DQD). (b) Scanning electron microscope im- 
age of the Si DQD before the top gate formation. The two side 
gates located next to side gate C are grounded, (c) Charge 
stability diagram of the Si DQD as a function of Vl and Vr 
at zero magnetic field, where Vds = — 2 mV, Vtg = 0.90 V, 
and Vc = —1.72 V. The white dotted lines are boundaries of 
the stable charge states. The charge numbers in the left and 
right QDs are Nl and N-r, respectively. 

mated to be 310 and 260 /ieV, respectively [22]. In confir- 
mation, AE can be approximated as AE = h 2 /87rm*A, 
where m* gives effective mass here, h is Planck's con- 
stant, and A is the area of the QD [23], with spin and 
valley degeneracies included. This equation determines 
AE to be between 260 and 380 fieV for our device geome- 
try [22] , in good agreement with the experimental estima- 
tion. We conclude that the QD is formed between the two 
constricted regions indicated by the ovals in Fig. [TJb). 

Current rectification in DQDs due to a PSB appears 
at a triple point with only one bias polarity [12 . We 
observed such current rectification with a negative bias 
voltage at triple point B in Fig.[TJc), as indicated by the 
trapezoid in Fig. [2|a), whereas no current rectification 
appeared with positive bias as shown in Fig. [2^b). In 
addition, the current rectification is lifted along the outer 
edge of the PSB regime indicated by the circle in Fig.|2|a) 
because of electron exchanges between the DQD and the 
right lead, comparable to PSB seen in GaAs DQDs [T2] . 

Since Si DQDs normally have doubly degenerate val- 
leys due to confinement in the direction perpendicular to 
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FIG. 2: (color online) (a) Triple point B shown in Fig. [TJc) 
with negative bias, where Vds = —2 mV, Vtg = 0.97 V, and 
Vc = —1.76 V. The PSB appears only for this polarity. Here 
e is the detuning axis, (b) The same triple point as in (a) 
under a positive bias (Vds = 2 mV). (c) Energy diagrams of a 
Si DQD at the circle marked in a (the left diagram) and at the 
blue cross marked in b (the right diagram), where the valley 
degeneracy is assumed to be lifted, (d) The same diagram as 
(c) without an assumption that lifting of the valley degeneracy 
is small. Intra-dot and inter-dot tunnelings between different 
valleys are assumed to be weak so that the PSB is not lifted. 

the Si surface, the valley degeneracy could lift a PSB. 
However, the fact that a PSB is observed indicates either 
a lifting of valley degeneracy or weak tunneling between 
valleys [24]. In the former case, once two spins occupy the 
(1,1) triplet state as shown in Fig.|2|c), the current flow 
is suppressed due to the PSB until relaxation from (1,1) 
triplet to (1, 1) singlet occurs. In the latter case, even if 
degenerate valleys exist as shown in Fig. [5|d), the PSB 
is not lifted because intra-dot and inter-dot tunnelings 
between valleys are weak. 

PSB features were observed at adjacent triple points, 
marked B, C, and D in Fig.JTJc). This is not expected for 
simple spin- \ PSB. Since the DQD has many electrons, 
spin- 1 ground states can exist, leading to scenarios for 
consecutive PSB [12 . Blockade where valley degeneracy 
plays a role can also lead to consecutive PSB-like fea- 
tures. Even when a spin doublet is formed in DQDs, the 
current flow could be suppressed because of weak tunnel- 
ing between valleys discussed above [22] . 

Figure j3^a) shows the leakage current in the PSB 
regime at triple point C in Fig. [TJc) as a function of 
magnetic field B applied normally to the DQD with a de- 
tuning, e, corresponding to the arrow shown in the inset. 
A strong current dip was observed at B = 0, whereas the 
current with opposite bias does not change as a function 
of magnetic fields [22]. Similar current dips have been 
observed for DQDs in In As nanowires p~4j [25] and car- 
bon nanotubes [15] and can be attributed to spin-orbit 
induced relaxation [26], which is suppressed at B = 



3 




FIG. 3: (color online) (a) Leakage current in the PSB regime 
as a function of magnetic field applied perpendicularly to the 
DQD and detuning, where V ds = -2 mV, Vtg = 0.97 V, 
and Vc = —1.99 V. Inset: Magnified plot of triple point C in 
Fig. [TJc), where the arrow corresponds to the detuning axis 
in the main figure, (b) Current along the dashed line in (a) 
denoted by the squares, and the fit to the data indicated by 
the blue line, (c) Values of Be extracted from the fit as a 
function of Vc- Large Vc corresponds to a large inter-dot 
tunnel coupling t. 

due to a Van Vleck cancellation p~4[ [27] . A Lorentzian 
line shape, J fit = i m ax{l - 8B^/9(B 2 + B 2 ,)} with char- 
acteristic width Be, is predicted theoretically [26 . The 
squares in Fig. |3|b) correspond to the absolute values of 
the leakage current in the PSB regime along the dashed 
line in Fig. |3^a). Fits to the Lorentzian form (the blue 
curve in FigTjSJb)) yield good agreement between theory 
and experiment. Furthermore, as the inter-dot tunneling 
between the two QDs is enhanced by changing Vc, the 
value of Bq extracted from the fit increases, as plotted in 
Fig. [3^c). This result is also consistent with the theory, 
which predicts Bq proportional to inter-dot tunnel cou- 
pling [26 . These results suggest that spin-orbit effects 
dominate spin relaxation in these devices. 

Another possible mechanism leading to a dip in current 
leakage around B = is spin- valley blockade with short- 
range disorder [28] , where the current dip as a function of 
magnetic-field-induced valley splitting is predicted. How- 
ever, we have no independent evidence that the required 
B-dependent valley splitting exists. The physics of the 
valley in Si DQDs deserves further experimental and the- 
oretical study. 

For some triple points, we observe a peak, rather than 
dip, in PSB leakage current on a larger field scale. As 
an example, the field dependence of the leakage current 
at triple point A in Fig. [TJc) is shown in Fig. |4ja). The 
arrow in the magnified plot of triple point A shown in 
Fig. [4^b) corresponds to the detuning axis in Fig. |4ja). 



FIG. 4: (color online) (a) Leakage current in the PSB regime 
as a function of magnetic field applied perpendicularly to the 
DQD and the detuning, whereVks = 2 mV, Vtg = 0.968 V, 
and Vc = —1.925 V. (b) Magnified plot of triple point A in 
Fig. mc), where the arrow corresponds to the detuning axis 
in (aj7 (c) Current along the dashed line in (a) denoted by 
the circles, and the fit to the data indicated by the blue line. 

Among the 15 triple points that show PSB [Fig. [ljc)], 
nine show a zero-field current dip and two show a peak. 
We also observed current peaks outside a current dip in 
some cases. 

In GaAs DQDs, zero-field peaks in leakage current 
were attributed to hyperfine-induced spin relaxation p~3j 
[29] . However, the contribution of the hyperfine interac- 
tion should be small in Si systems, because the dominant 
28 Si atoms have zero nuclear spin. Using 4.7% natural 
abundance of 29 Si and lithographic device dimensions [22] 
gives an expected number N of nuclear spins in a Si DQD 
to be 2-3 x 10 4 , corresponding to a fluctuating Over- 
hauser field magnitude B nuc = \A\/ gii^yfN ~ 10-15 //T, 
where the hyperfine coupling constant \A\ ~ 0.2 /ieV from 
NMR measurements [31 and g ~ 2 for electrons in Si. 
Since the peak width in Fig. |4^c) is larger than B nuc by 
a factor of 10 4 , the mechanism of the current peaks at 
B = is not explained by hyperfine interaction. 

Similar peaks were also seen in Si DQD in Ref. [19], 
where the peak is well described by spin-flip cotunneling 
[32] . When k^T > t (&b is Boltzmann's constant and t is 
the inter-dot tunnel coupling), the spin-flip cotunneling 
current is given by I cot = Aecgji^B /3sinh(g/iB^ /k&T) 
withc = /i[(r R /(A-e)) 2 + (r L /(A + e-2/7 / -2eF ds )) 2 ]/^ 
where T L (r) is the coupling of the lead to the left (right) 
dot, A is the depth of the two-electron level [33], and U' 
is inter-dot charging energy. Since we observed clear reso- 
nant tunneling peaks, T L (r) is larger than t [34 . In addi- 
tion, if r L (R) > t > ksT ~ 21 /ieV, the current would be 
much larger than the observed current shown in Fig.|4^b). 
As a result, k^T > t so that I cot can be used to fit the cur- 
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rent peak. The blue curve in Fig. |4jc) is 7 cot , which has 
a good agreement with the data by using T ~ 250 mK, 
yielding g ~ 2.3 and c ~ 54 kHz//ieV. Since the cur- 
rent does not vary much along the base of the triangle in 
Fig. |4jb), we assume Tl ~ Tr = T. By using expression 
of c with A ~ 1 meV, e ~ meV, U' ~ 1 meV, and 
e Vd s ~ 2 meV estimated from the bias triangle shown in 
Fig. [4jb), we extracted T ~ 26 /ieV. Furthermore, £ can 
be extracted to be about 0.3 jaeV from the unblocked res- 
onant tunneling peak current (~ 0.6 pA) with Eq. (15) in 
Ref. [21] . These values are similar with those in Ref . [T9] 
and in an experimentally reasonable range so that the 
spin-flip cotunneling processes are most likely the mech- 
anism of the peak. It should be noted that, as for the 
dip in Fig. 3, spin- valley blockade with disorder could 
also explain the peak, but again we have at present no 
evidence of the required field-dependent valley splitting 
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